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<— | ! Abstract 

The QED vacuum polarization function is calculated to C(a 2 ) (next-to- leading order) accuracy in 
k> , the threshold regime by using the concept of effective field theories to resum diagrams with the 

instantaneous Coulomb exchange of longitudinally polarized photons. It is shown that the 0(a 2 ) 
contributions are of order a 2 in size rather than a 2 /it 2 . The vacuum polarization contributions to 
the 0(a 6 ) hyperfine splitting of the positronium ground state are recalculated and differences with 
an older calculation are pointed out. The results are used to determine 0(Cpa 2 ) (next-to-next- 
to-leading order) Darwin corrections to heavy quark-antiquark bound state I — wave functions 
at the origin and to the heavy quark-antiquark production cross section in e + e~ collisions in the 
threshold region. The absolute value of the corrections amounts to 10% — 20% and 17% — 34% in 
the modulus squared of the ground state wave functions at the origin for the bb and cc systems, 
respectively. In the case of the ti production cross section in the threshold region the absolute value 
of the corrections is between 2% and 6% around the 15* peak and between 1% and 2% for higher 
energies. A critical comment on recent QCD sum rule calculations for the T system is made. 



1 Introduction 

In recent years many sophisticated methods have been developed to calculate higher order ("multi- 
loop") QCD radiative corrections for high energy quantities for which it is believed that an expansion 
in terms of Feynman diagrams with a certain number of loops represents an excellent approximation 
to the predictions of quantum chromo dynamics. Notable examples are the hadronic cross section in 
e + e~ collisions at LEP energies or the (photonic) vacuum polarization function. In the high energy 
limit, where the quarks can be treated as massless, these quantities have been calculated up to three 
loops (l|, H], m ||]. However, future experiments (NLC, B-factory and r-charm factory) will test the 
vacuum polarization function and the hadronic cross section also in the kinematic regime close to 
heavy quark-antiquark thresholds, where bound state effects become important. The threshold regime 
is characterized by the relation 




\0\ <a s , p= Wl-4^-, (1) 

V q z + te 

where Mq is the heavy quark mass and vT denotes the cm. energy. In the process of heavy 
quark-antiquark production above threshold, q 2 > AMq, (3 is equal to the velocity of the quarks in 
the cm. frame. We therefore call (3 "velocity" in the remainder of this work, even if q 2 < &Mq. 
In the threshold regime the accuracy of theoretical predictions to the hadronic cross section and to 
the vacuum polarization function is much poorer than for high energies. Aside from definitely non- 
perturbative effects (in the sense of "not calculable from first principles in QCD"), the breakdown of 
the perturbative expansion in the number of loops makes any theoretical description in the threshold 
region difficult. This breakdown of the perturbation series is indicated by power (1//3) or logarithmic 
(ln/3) divergences in the velocity which blow up if evaluated very close to the threshold point. Some of 
these divergences (e.g. the a n / ' (5 n , n > 1, Coulomb singularities in the Dirac form factor iq describing 
the electromagnetic vertex) can be treated by using well-known results from non-relativistic quantum 
mechanics, but a systematic way to calculate higher-order corrections in the threshold regime seems 
to be far from obvious, at least from the point of view of covariant perturbation theory in the number 
of loops. This type of perturbation theory will be referred to as "conventional perturbation theory" 
from now on in this work. 

On the other hand, there are many examples of heavy quark-antiquark bound state properties 
where the complete knowledge of higher-order corrections would be extremely valuable. Most of the 
present analyses (see e.g. ||) are based on leading and next-to-leading order calculations. Here, higher- 
order corrections could significantly increase the precision of present theoretical calculations, but could 
also serve as an instrument to test how trustworthy certain theoretical predictions are and to estimate 
the size of theoretical uncertainties. Further, they could contribute toward a better understanding of 
the role of non-perturbative effects (in the sense mentioned above) in apparent discrepancies between 
the determination of the size of the strong coupling from the T(IS') decay rates || and QCD sum rule 
calculations for the T system [pJF] on the one hand, and from the LEP experiments on the other. 

The framework in which bound state properties and also dynamical quantities in the threshold 
regime can be calculated in a systematic way to arbitrary precision is non-relativistic quantum chro- 

1 During completion of this paper we became aware of a new publication, where QCD sum rules for the T system are 
used to determine the strong coupling and the bottom quark mass [H]. We will give a brief comment on this publication 
and on [|7| at the end of Section ||5[ 



modynamics (NRQCD) || which is based on the concept of effective field theories. In the kinematic 
regime where bound states occur and slightly above the threshold, NRQCD is superior to conventional 
perturbation theory in QCD and (at least from the practical point of view) also to the Bethe-Salpeter 
approach, because it allows for an easy and transparent separation of long- and short-distance physics 
contributions. This is much more difficult and cumbersome with the former two methods. However, 
we would like to emphasize that all methods lead to the same results. As an effective field theory, 
NRQCD needs the input from short-distance QCD in order to produce viable predictions in accor- 
dance with quantum chromodynamics. This adjustment of NRQCD to QCD is called the matching 
procedure and generally requires multi-loop calculations in the framework of conventional perturbation 
theory at the level of the intended accuracy. 

In this work we demonstrate the efficient use of the concept of effective field theories to cal- 
culate the QED vacuum polarization function in the threshold region to 0(a 2 ) accuracy. In order 
to convince the reader of the simplicity of the approach we use our result to recalculate the vacuum 
polarization contributions to the 0(ofi) hyperfine splitting of the positronium ground state energy 
level without referring back to the Bethe-Salpeter equation. Differences between our result and an 



older calculation ||, 1C], are pointed out. We analyse the vacuum polarization function at the bound 
state energies and above threshold and, in particular, concentrate on the size of the 0(a 2 ) corrections. 
It is shown that the size of the 0(a 2 ) corrections in the threshold regime is of order a 2 rather than 
a 2 /it 2 which is a consequence of their long-distance origin. In a second step our results for the QED 
vacuum polarization function are applied to calculate 0(C F a 2 ) (next-to-next-to-leading order) Dar- 
win corrections to the heavy quark-antiquark I = bound state wave functions at the origin and to 
the cross section of heavy quark-antiquark production in e + e~ annihilation (via a virtual photon) in 
the threshold region. The corresponding unperturbed quantities are the solutions of the Schrddinger 
equation for a stable quark-antiquark pair with a Coulomb-like QCD potential, VQ CD (r) = —Cpa s jr. 
It is demonstrated that the size of the 0(C F a 2 ) Darwin corrections is also of order a 2 rather than 
a 2 /7r 2 . We present simple physical arguments that the scale of the strong coupling governing the 
0(C F a 2 ) Darwin corrections is of order CfMqcx s and we analyze the size of the corrections for the 
tt, bo and cc systems assuming that the size of the Darwin corrections can be taken as an order of 
magnitude estimate for all (yet unknown) 0(a 2 ) corrections. The sign of the latter corrections and 
their actual numerical values can, of course, only be determined by an explicit calculation of all 0(a 2 ) 
contributions. 

At this point we want to emphasize that our approach does not depend on any model-like 
assumptions, but represents a first principles QCD calculation. The only assumptions (for heavy 
quarks) are that (i) the instantaneous (i.e. uncrossed) Coulomb-like exchange of longitudinal gluons 
(in Coulomb gauge) between the heavy quarks leads to the dominant contributions in the threshold 
regime and is the main reason for heavy quark-antiquark bound state formation and that (ii) all further 
interactions can be treated as a perturbation. We believe that the actual size of the 0(a 2 ) corrections 
can then serve as an important a posteriori justification or falsification of these assumptions for the 
different heavy quark-antiquark systems. Finally, we address the question whether bound state effects 
can lead to large corrections to the vacuum polarization function in kinematical regions far from the 
actual threshold regime. We come to the conclusion that such corrections do not exist. 

The program for this work is organized as follows: In Section [2] the calculation of the QED 
vacuum polarization function to 0(a 2 ) accuracy in the threshold region is presented. We define a 
renormalized version of the Coulomb Green function for zero distances, which allows for application 



of (textbook quantum mechanics) time-independent perturbation theory to determine higher-order 
corrections to wave functions and energy levels. For completeness we also give an expression for the 
QED vacuum polarization function valid for all energies with 0(a 2 ) accuracy. In Section the QED 
vacuum polarization function in the threshold region is analysed with special emphasis on the size of 
the 0(a 2 ) corrections, and the 0(a e ) vacuum polarization contributions to the positronium ground 
state hyperfine splitting are calculated. Section || is devoted to the determination and analysis of 
the 0(Cpa 2 ) Darwin corrections to the bound state wave functions at the origin and the production 
cross section in the threshold regime for the different heavy quark- antiquark systems. In Section || we 
comment on the existence of threshold effects far from threshold and Section y contains a summary. 



2 Determination of the QED Vacuum Polarization Function 
in the Threshold Region 

We consider the QED vacuum polarization function II defined through the one-particle-irreducible 
current-current correlator 

( g V - gV)n(« 2 ) = -i fd 4 xe^ x {o\Tf(x)f(o)\o), (2) 



where j^(x) = ie^>{x)^^>(x) denotes the electromagnetic current. V? represents the Dirac field of the 
electron with charge e. According to the standard subtraction procedure, II vanishes for q 2 = 0. It 
has been shown in |11] that in the kinematical region close to the e + e~ threshold point q 2 = AM 2 , 



where M denotes the electron mass, the current-current correlator II is directly related to the Green 
function G° E {x,x') of the positronium Schrodinger equation 



M x \x\ 



G%{x,x') =5^(x-y), (3) 



where E denotes the energy relative to the threshold point, E = y/q 2 — 2M, and a is the fine structure 
constant. Explicit analytic expressions for the Green function have been calculated in a number of 
classical papers [[E|. The proper relation between the vacuum polarization function and the Green 
function in the threshold region reads [11] 

nr»tf) = ^4(o,o). (4) 

Because we are only interested in 0(a 2 ) accuracy in the threshold region we effectively can replace 
the factor 1/q 2 in eq. (f§) by 1 /AM 2 . 

For illustration, let us now examine the one- and two-loop contributions to the vacuum polar- 
ization function and the expression for the vacuum polarization function from non-relativistic quantum 
mechanics according to eq.(||). The one- and two-loop contributions to II have been known for quite 
a long time for all energy and mass assignments fll3| , [14| , |lq| . Far from the threshold point those loop 
results provide an excellent approximation to the QED vacuum polarization function at the 0(a 2 ) 
accuracy level, 

iCoV) = (^KV) + (£) a n<*V) + o(Q 3 ) . (5) 



In the kinematic domain where a <C \(3\ <C 1, the expansion in terms of the number of loops is still an 
adequate approximation, and we are allowed to expand the coefficients in eq. (J5|) for small velocities, 
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+ 0(a 2 (3), (7) 



where = a/1 -4M 2 /(g 2 + ie) and (3 = 1.202056903 .. .. In eq. © the £>(a/3) contribution is also 
displayed, allowing for a check of the normalization during the matching procedure. Whereas the one- 
loop contribution, eq. (J6(), can be evaluated for — ► 0, indicating that the O{o) contribution of the 
vacuum polarization function is of pure short-distance origin, the two-loop expression, eq. (|7|), diverges 
logarithmically for vanishing (3. This shows that beyond 0{a) accuracy conventional perturbation 
theory is inadequate for |/3| < a due to long-distance effects which cannot be calculated in terms of 
a finite number of loop diagrams. These long-distance effects can, on the other hand, be described 
adequately by the vacuum polarization function calculated in the framework of quantum mechanics, 
which per constructionem is valid in the non-relativistic limit. However, the vacuum polarization 
function as defined in eq. @ gives a divergent result ,p 
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where 7 is the Euler constant and ^ represents the digamma function, 
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i=l 



0.5772156649. 



i&(z) = — lnT(z). 

The divergences in eq. (J8[) can be easily understood if non-relativistic quantum mechanics is considered 
as a low-energy effective theory which is capable of describing long-distance physics close to the 
threshold (characterized by momenta below the scale of the electron mass) but does not know per se 
any short-distance effects coming from momenta beyond the scale of the electron mass. This lack of 
information is indicated in eq. (JSJ) by short distance (UV) divergences and has to be cured by matching 
non-relativistic quantum mechanics to QED. The result of this matching procedure is called "non- 
relativistic quantum electrodynamics" (NRQED) |l6|]. In this light we have to regard relations (||) 
and (0) as unrenormalized, which we have indicated by using the superscripts 0. 

In the common approach the Lagrangian of NRQED is obtained by introducing higher dimen- 
sional operators in accordance with the underlying symmetries of the theory and by matching them to 



only. 



In eq. (H) we can identify /3 = yl — 4M 2 /(g 2 + ie) with *J{E + ie)/M because we are interested in 0(a 2 ) accuracy 



predictions within conventional perturbation theory in QED in the kinematical regime a <C |/?| "C 1 
where both NRQED and QED are valid and must give the same results^. In general this leads to 
divergent renormalization constants multiplying the NRQED operators which then cancel the diver- 
gences in eq. (||) and add the correct finite short-distance contributions. In our case, however, the 
explicit determination of the renormalization constants is not necessary because we can match the 
vacuum polarization function obtained from unrenormalized NRQED directly to the one- and two- 
loop expressions from QED, eqs. (J6|) and (0). This "direct matching" method has the advantage that 
the regularization of the UV divergences in eq. (|8|) can be performed in a quite sloppy way, but has the 
disadvantage that it is of no value to determine other quantities than the vacuum polarization function 
itself. To arrive at the vacuum polarization function in the threshold region we just have to replace 
the /3-independent and divergent contributions of the 0(a) and 0(a 2 ) coefficients of expression 
an expansion for small a 
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by the corresponding 0- independent and finite contributions in eqs. ©, (0)- The correct normaliza- 
tion between the contributions coming from the loop calculations and from non-relativistic quantum 
mechanics can be checked explicitly by observing that the coefficients in front of the /3-dependent 
0(a) and 0(a 2 ) contributions in eqs. (|6|), (^) and ([]) are identical. At this point we would like to 
note that the function A contains only contributions of order a 3 and higher for an expansion in small 
a. We shall return to this point later. The final result for the vacuum polarization function valid to 
0(a 2 ) accuracy in the threshold region then reads 
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It is an interesting fact that this result can be obtained directly from the one- and two-loop results, 
eqs. @ and ([?]), by the replacement 



= \n(-i f3) - ^ A(a, P) . 



(12) 



As mentioned earlier, the function A is of order a 3 for a ^C \/3\. In the language of Feynman diagrams 
A arises from diagrams with instantaneous Coulomb exchange of two and more longitudinal photons 



For a <C |/3| <C 1 conventional perturbation theory in QED is valid because a represents the smallest parameter, 
whereas NRQED is valid because \/3\ is much smaller than the speed of light. 



(in Coulomb gauge) between the electron-positron pair. However, if \/3\ is smaller than a {(3 being 
real), then A is of order a 2 , 



A(a,/3) 



\/3\<a a 



2 r 



In [i 



2/3 



n 



7 



+ 0{a(3) 



(13) 



At this point it is illustrative to examine the limits a <C /3 and /3 <C a for the function H, defined in 
eq. fll2|) , for real and positive values of j3: 



H(a,(3) a ^ \n(-if3) + 0Q, 
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It is evident that the function H interpolates between a In /3-behaviour in the region where conventional 
perturbation theory is valid and a constant with a logarithm of a for /3/a — > 0. This leads to a finite 
value for II T ^ at the threshold point. As we will see in the next section, H T f^, has singularities 
at the positronium energy levels indicating that the breakdown of conventional perturbation theory 



is directly related to the formation of bound states of the virtual e + e~ pair [11|. 

Based on result (fiA\ ) we are now able to define a renormalized expression for the zero-distance 
Green function by inverting relation (Eh 
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As we will show later, this renormalized zero-distance Green function can be used for the calculation 
of higher-order corrections in time-independent perturbation theory. 

For completeness we also present the QED vacuum polarization function with 0(a 2 ) accuracy 
for all energies, 
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where 



1-/3 
1 + /9 



and Li2, Li3 denote the di- and trilogarithms [17]. The reader should note that IIqeS vanishes at 



q 2 = and is an analytic function in q 2 except at poles and branch cuts, and satisfies the dispersion 
relation 
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The explicit form of ImilqED in the threshold region will be presented in Section |3|. For the use and 
interpretation of formula ( |T7| ) see also Section 0. 



3 Examination of the Vacuum Polarization Function 
in the Threshold Region 

In this section we analyse the properties of the vacuum polarization function in the threshold region 
above and below the threshold point q 2 = AM 2 . Compared to an older work on the same subject [11] 



we are not so much interested in general properties of perturbation theory in the presence of bound 
state formation, but in the explicit form and behaviour of Hp^" . In particular, we focus on the size 
of the 0{a 2 ) contributions. We also would like to mention that the vacuum polarization function 
has been studied in a similar way in ||. In the latter publication, however, a different definition for 
the vacuum polarization is employed, only the positronium ground state energy is considered and a 
contribution in the one-loop result is missing. We will come back to this point later. Comparing 
the methods used in || with the effective field theoretical approach employed in this work makes the 
elegance of the latter technique obvious. 

We start in the kinematic region above threshold where a <C (3 C 1 Here, as mentioned in 
the previous section, the one- and two-loop results, eqs. ©-(0), are reliable. This is consistent with 
the fact that the function A contains only contributions of order a 3 and higher, 
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Thus, for practical applications the contributions of function A can be neglected. (See also the dis- 
cussion in Section |5|.) One might think that for a ~ (3 the one- and two- loop expressions should still 
represent an appropriate 0(a 2 ) prediction, because the radius of convergence of the series on the r.h.s. 
of eq. (|l9|) is \/3\ = a/2 |l8|]. However, as illustrated in Fig. [I], for a ~ (3 the contributions coming from 
function A are already of order a 2 /n 2 and thus have to be included if 0(a 2 ) accuracy is intended. 
For even smaller velocities, of course, the contributions from A are essential because they cancel the 
divergent In/? term from the two-loop expression n*- 2 ^, see eqs. (]?]) and (|l3|). Therefore the value of 
^Thr & t t ne threshold point is finite and reads^ (a = 1/137) 

The plus sign in the argument of IT T ^" indicates that the expression on the r.h.s. of eq. (|20|) represents only a 
right-sided limit on the real q 2 -axis. 
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Figure 1: The 0(a 2 ) corrections to the vacuum polarization function in the threshold region with and 
without the contributions contained in function A, eq.dTH), in the kinematic region < (3 < 2a above 
the threshold. The solid line denotes Re [IT 2 ) + it 2 /a 2 A], the dashed line RelT 2 ), the dashed-dotted 
line Im[IT 2 ) + ir 2 /a 2 A] and the dotted line ImlT 2 ). The value of the fine structure constant is taken 
as a = 1/137. IT 2 ) represents the two-loop contribution to the vacuum polarization function and is 
displayed for in eq. ((J). 
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(20) 



It is evident from eq. (pc|) and Fig. [T] that the size of the 0(a 2 ) corrections in the threshold region is of 

order a 2 rather than a 2 /ir 2 , whereas the O{o) contribution is of order a/ir. This can be understood 

from the fact that the 0{a) contribution in IXt^ comes entirely from the one-loop result IT 1 ), 

eq. (||), and therefore originates from momenta beyond the scale of the electron mass. High momenta 

contributions are expected to be of order a/V if no "large logarithms" occur^j. The large 0(a 2 ) 

contributions to II T ^ , on the other hand, arise from the interplay of the logarithm of the velocity 

in IT 2 ), eq. (0), and the contributions from the instantaneous Coulomb-exchange of two and more 

longitudinal photons between the virtual electron-positron pair. For small velocities the latter effects 

generate a logarithm of the velocity with an opposite sign, which cancels the logarithm in IT 2 ) . We 

5 For comparison the reader might consider the well-known one- and two-loop contributions to the anomalous magnetic 
moment of the electron Q [jlj, g e - 2 = (f ) - 0.66 (f ) 2 + C((f ) 3 ). Here, long-distance effects from the e+e" threshold 
play no role. Therefore g e — 2 can be regarded as a typical short-distance quantity with no "large logarithms" . 



therefore conclude that the large 0(a 2 ) contributions are of long-distance origin. This is particularly 
obvious for the In a term which could never be generated at short distances. 

The situation for a <C \j3\ <C 1 below threshold is similar to the one above threshold. Here, 
the one- and two-loop contributions from conventional perturbation theory, eqs. (|5J)-(0), provide a 
viable prediction, because the contributions from the function A are of order a 3 and higher. They are 
beyond the intended accuracy and can be neglected. (See also the discussion in Section |5[) On the 
other hand, it is obvious that the one- and two-loop results are not sufficient for energies close to the 
positronium bound state energies, 
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because the vacuum polarization function is expected to have poles at those energy values. Therefore 
the full expression for n T ^" , eq. (|TT1), must be employed. It is straightforward to check that n T ^" 



indeed has poles at the positronium energy levels [11] leading to the following Laurent expansion at 
the bound state energies E n = — . % , (n = 1, 2, 3, . . .), 
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For completeness we also present the corresponding Laurent expansion for the renormalized zero- 



distance Green function based on definition (16), 

|2 A T\f2 



lim G|(0,0) 



l*n(0)|' 



E n — E 



ie 



4 M 1 
9^" 



+ 



M 2 a 
2ir 



+ 0(E n - E) 



(24) 



As expected, the residues at the bound state energies are equal to the moduli squared of the normalized 
I = Coulomb wave functions at the origin, 
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In eqs. (|22D-(|24|) we have also displayed the constant terms of the Laurent expansion. These constants 
are relevant for higher-order corrections to the positronium energy levels and to the wave functions at 
the origin. The size of the 0(a 2 ) corrections in these constant terms is (similar to the 0(a 2 ) contribu- 
tions above threshold) of order a 2 rather than a 2 /it 2 indicating again the long-distance character of 
the 0(a 2 ) corrections. In Table |] we have displayed the numerical values of a n for the radial quantum 
numbers n = 1, 2, 3, 4, 5. It is an interesting fact that the n — ► oo limit of a n exists 
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and coincides with 0(a ) contributions of ReII T ^ r '{q — > AM +), eq. (PQ). The numerical value for 
linin^oo a n is presented in Table |l[ 

To illustrate the importance of the constants a n in time-independent perturbation theory 
(TIPT), we recalculate the O(o 6 ) "vacuum polarization" contributions to the ground state triplet- 
singlet hyperfine splitting (HFS) of the positronium, which were, to our knowledge, considered for 
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n 


1 


2 


3 


4 


5 


oo 


a n 


2.89 


2.48 


2.35 


2.29 


2.25 


2.10 



Table 1: The numerical value for the constants a n for the radial quantum numbers n = 1,2,3,4,5 
and for n — ► oo with a = 1/137. 



the first time in ||. The vacuum polarization contributions to the HFS in the energy levels of the 
positronium system arise from the effect that the bound triplet ( 3 <Si, J = 1 ) e + e~ pair can 
annihilate into a virtual photon for a time period of order 1/iVf , whereas the singlet ( 5o, J = 0~ + ) 
cannot. If the virtual photon energy is approximated by \Jlp = 2M, this annihilation process leads 
to a (5-function kernel in the coordinate-space representation (corresponding to a constant kernel in 
momentum space) with the form, 

2an 



HAnnix) 



M 2 



■<5®(x) 



(27) 



r O( Q 2) 



This kernel can now be used in TIPT. Taking into account that 11^" ; contains 0(a) as well as 
0{a 2 ) contributions we have to apply second- and third-order TIPT to obtain all relevant 0(a e ) 
contributions to the HFS. The formal result for the 0(ce 6 ) energy shift for the triplet states with 
radial quantum numbers n and with I = due to Hauu reads 



rpa" 
0Cj Ann,n 



\ n I H Ar , 



E n — El 



H Ann | n ) 



l^n 



+ 



t t< 



n | H Ar . 



m) (to I 



HAr. 



k)(k 



HAr. 



n 



0{cfi) 



(28) 



m^n k^n 

where | i), i = I, to, n, k, represent normalized (bound state and free scattering) eigenfunctions to the 
positronium Schrodinger equation with the eigenvalues £"j. The symbol {}o(a (i ) indicates that only 
0(a 6 ) contributions are taken into account. It is evident from the form of HAnn(%) that only the 
zero-distance Green function is relevant for 5E^ nnn , 



l+n 



Ei — E„ 



0) 



l^n 



l^(Q)l s 

Ei - E n 



lim 

E—*E n 



G° E (0,0) 



l*n(0)| 



E n — E 



l€ 



(29) 



However, relation ( p9[ ) still contains divergences (see eq.(||)). As we have pointed out in Section g, these 
divergences indicate that non-relativistic quantum mechanics is not capable to describe physics if the 
relative distance of the electron-positron pair is smaller than the inverse electron mass. Therefore, we 
have to replace G° E (0, 0) in relation (29) by its renormalized version G E (0, 0), eq. ( |I6| ) (using eq. (11)), 



which describes short-distance physics properly. The final expression for 5E'^ nn n then reads 
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Ma 6 f 1 , 1/1 , ^-! 1 \ 1/13 21, \ 11 I 



- In a- - - + lnn- V - + — - — + — (3 + - In2 L (30) 



4n 3 12 ' 2\n frf i J 4tt 2 V81 2*/ 32 4 

Taking also into account that the virtual photon energy is smaller than 2M by the amount of the 
binding energy E n = — 4 ^> we have to replace M 2 in eq. (^?j) by (M + E n /2) 2 and therefore get 
another C(a 6 ) contribution, 

4^ ' Ion 

where 6E\ nnn represents the C(a 4 ) energy shift due to Hauu- For the ground state (n = 1) the 
complete C(a 6 ) vacuum polarization contribution to the HFS then reads 

«BSL. 1 + «tf*. " ^T ! {5 ll " , + 4^(I + l^)-i + I 1112 }- (32> 

Our result differs from the one presented in ||, [l(| by the amount Ma 6 /8. Half of the discrepancy comes 
from the fact that in M the binding energy contribution 5Eg ind 1 was not taken into account, whereas 
the other half originates from a missing 0(a 2 ) contribution in the one-loop vacuum polarization^]. 

Before we turn to applications of our results in the context of QCD, we do not want to leave 
unmentioned that the leading contributions to the normalized cross section for production of a heavy 
lepton-antilepton pair (with lepton mass M) in e + e~ collisions (via a virtual photon) in the threshold 
region can be recovered from IlyL. by means of the optical theorem |2C, J2jJ , 



& - _ „(.-«+^L+L-) _ l ImlgV) _ gL Im0g , 



= J? g l^(0)l 2 HB - E„) + S(E) - 1 _ J^_ V) , (33) 

where a p t represents the point cross section and only final-state interactions are taken into account. 



4 Darwin Corrections in QCD 

In the previous section we have shown that the size of the 0(a 2 ) corrections to the QED vacuum 
polarization function in the threshold region is of order a 2 rather than a 2 /it 2 . Although this fact is 
important for precision tests of QEDFl it does not lead to theoretical concerns about the convergence 

6 This can be easily seen by comparing eq. (33) of |9] with eq. (pi) in this work for 

2 /„ A/f Ma 2 2 a ■ a , ml 3\ 

This shows that at bound state energies the one-loop contribution to the vacuum polarization function also contains 
terms of order a 2 . 

7 As far as tests of QED in the t + t" system in the threshold region are concerned the present experiments do not 
even reach the O(a) (next-to-leading order) accuracy level. This can be easily seen from the fact that the complete 
threshold region for the t + t~ system, |/3| <J a <=>■ \\/<p — 2m T \ < m T a 2 = 0.1 MeV still lies within the limits on the 
tau mass itself, m T — 1777.001q'27 MeV |22|. Thus only experiments on electron and muon systems can be regarded as 
precision tests of QED in the threshold regime. 
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of the perturbative series because of the smallness of the fine structure constant a and because QED 
is not asymptotically free. 

In the framework of QCD, however, the situation is completely different: the coupling is much 
larger and even becomes of order one for scales much lower than 1 GeV. Therefore, the fact that 
the size of the 0{a 2 s ) (next-to-next-to-leading order) corrections in the threshold region might be of 
order a 2 rather than ot 2 /ir 2 is an extremely important theoretical issue because this would lead to 
corrections of order 1% — 25% rather than 0.1% — 2.5% for a s = 0.1 — 0.5. Here, two natural questions 
arise: what scale should be used in the strong coupling, and for which heavy quark-antiquark systems 
the 0{a 2 ) corrections represent contributions to the asymptotic perturbation series in the convergent 
regime. These questions will be addressed in the following section. 

To be more specific we will calculate the 0{Cpa 2 ) Darwin corrections to the (5-wave, I = 0) 
wave functions of a bound heavy quark-antiquark pair at the origin and to the heavy quark-antiquark 
pair production cross section in e + e~ annihilation (via a virtual photon) in the threshold region. A 
presentation of all 0{Cpa 2 ) corrections including all kinematic and relativistic effects will given in a 
subsequent publication. The corresponding uncorrected quantities are the well-known exact solutions 
to a pure Coulomb-like non-relativistic quark-antiquark system described by a Schrodinger equation 
with the QCD-potential Vqcd{t) = -C F &■, where C F = (N 2 - 1)/2N C = 4/3. 

The Darwin interaction is generated in the non-relativistic expansion of the Dirac equation. 
In the coordinate-space representation it is proportional to a 5-function and reads[] 

H Dar (x) = ^f^S^(x). (34) 

A practical application for the corrections to the bound state wave functions at the origin is the 
leptonic decay rate of the J/^f and the T(15) and (maybe) of the first few excited states of the 
T family, whereas the corrections to the cross section would be relevant for tt production at NLC. 
We explicitly mention those applications in this context because it is believed that for them non- 
perturbative (in the sense "not calculable from first principles in QCD") effects are either well under 



control or even negligible [23, 20|. But, of course, these corrections can be applied to other heavy 
quark-antiquark systems as well, at least in order to check their size. At this point we want to 
emphasize that we do not intend to present a thorough phenomenological analysis in this work. The 
primary aim is to use the <D{Cpa 2 ) Darwin corrections to illustrate the typical size of the complete 
(and yet unknown) 0(a 2 ) corrections for the tt, bb and cc systems. Their actual numerical value and 
even their sign cannot, of course, be predicted at the present stage. 

To keep our analysis transparent we ignore all 0(a s ) corrections, the effects from the running 
of the strong coupling and also non-perturbative contributions like the gluon condensate. The latter 
effects are well known and have been treated in a large number of earlier publications. We further 
neglect the width of the quarks and treat them as stable particles for the most part in the following 
analysis. From the technical point of view the calculations of the 0(Cpa 2 ) Darwin corrections are 
identical to the corresponding QED calculations, which means that we use time-independent pertur- 
bation theory. However, we have to take care about the correct implementation of the number of 
colors, N c = 3, and the group theoretical factor Cf- In the following the superscript "QCD" indicates 
that the corresponding quantity is obtained from the QED expression by the replacement a — ► CpOi s . 

8 Compared to the Darwin interaction known from the hydrogen atom the expression on the r.h.s. of eq. (B4j) is a 
factor of two larger because both quark-antiquark-gluon vertices contribute. 
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It is then straightforward to determine the 0(Cpa 2 ) Darwin corrections to the modulus squared of 
the I = bound state wave functions at the origin, (n = 1,2,3,.. .), 



swrmh* 



-2|^ CD (0)H^^ lira 



m CD m 2 C 2 F aia 



M Q E->E. 
-! / •-! „ 2 -QCD 
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QCD 



E-ie 



(35) 



where the symbol {}o(c 2 a 2 ) indicates that only 0(C F a^) corrections are taken into account^. The 
calculation of the 0{Cpa 2 s ) Darwin corrections to the quark-antiquark cross section in the threshold 



region is more involved. Here, we apply the optical theorem, eq. (33), to the corrections of the 
zero-distance Green function themselves, 
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(36) 



The 0(Cpaf) Darwin corrections to the cross section above threshold then read 
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(37) 



where Rp^ r represents the "Sommerfeld factor" (sometimes also called the "Fermi factor"} 
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(38) 



Below threshold we have to determine the corrections to the residues of G E D (0, 0) at the bound 
state energies, where as shown above the corresponding bound state poles have to be subtracted. This 
calculation is straightforward and leads to the corrections to the I = bound state wave functions 
at the origin presented in eq. (^). It is an interesting fact that eq. ([??]) allows for the calculation of 
the shifts of the QQ bound state energies due to the Darwin interaction. To show this we rewrite the 
sum of the Sommerfeld factor, eq. (38), and the contribution involving the digamma function of the 
0(Cpa 2 ) Darwin corrections above threshold, see eqs. fllQ|), ( JTT| ) and ([37]), as 
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(39) 



9 Eq. ( |3q ) also generates 0(Cfo.s) corrections which differ from the well-known 0{Cfo/. s ) corrections generated by 
the (1 — ACfols/it) correction factor |24]. Adding up all the 0(Cfol s ) corrections and the corresponding renormalization 
constants will of course yield the correct result. The same remark holds for the result for the cross section above threshold, 

eq. 



©• 
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It can be easily checked that the function r( F £* a + 1 — i %g s ) develops poles at the energies^ 

£Qcd _ ^qcd + 6E ^ r t (n = 1, 2, 3, . . .) , (40) 

where the 5E® C p ar represent the energy shift of the I = Coulomb energy levels with the radial 
quantum number n generated by the Darwin interaction, 



■•'-" lf / ,,'-" ]i i rr. i ,.'" ■••■• - hi/'.'' "/ml 2 C Fa s ir 

M 2 Q 



SETDar = (">** I H Dar | n QCD ) = |^ CD (0)|' 



M Q C*a\ (41) 



!n 3 



At this point we also want to emphasize that the ln(C^a s ) and digamma contributions occurring 
in eqs. ( p5| ) and (K7l) are not related to the running of the strong coupling and therefore cannot be 
resummed by any known type of renormalization group equation in the sense of a leading logarithmic 
resummation. These logarithmic terms arise because two scales are relevant in the threshold region, the 



heavy quark mass Mq and the relative momentum of the quark- antiquark pair oc Cf ol s Mq [25|. The 
ln(Cj7 a s ) contributions induced by the running of the strong coupling have been discussed in [26, |23|| . 
Before we turn to the discussion on the size of the 0(C F a 2 s ) Darwin corrections we have to 
address the question of which scale one should use in the strong coupling. Strictly speaking, a final 
answer to this problem would require an O(a^) analysis, which is beyond the scope of this work. 
However, one can find simple arguments that the scale in the strong couplings of expressions (|35|) and 
(|37|) should be of the order CfO-sMq, which will be called "the soft scale" in the remainder of this 
work. We would like to remind the reader that the scale of the strong coupling in the unperturbed 
(pure Coulomb) quantities |^ / ,^ CD (0)| 2 and Rf£ r ls °^ the or der of the soft scale. This is obvious 
for the wave functions of the ground state and the first few excited states and for the cross section 
in the kinematic region (3 ~ Cfoc s because they describe bound quark-antiquark pairs with relative 
momentum of order Cf o s Mq. But it is also true for highly excited states (n 3> 1) and the cross 



section right at the threshold due to "saturation" effects |26|, [23[, i.e. the scale of the strong coupling 
is of order Cfo s Mq although the kinematical relative momentum of the quark pair vanishes^]. To 
understand that the scale of the 0(C F a 2 ) Darwin corrections should also be of order of the soft scale, 
let us have a closer look at the origin of the strong couplings governing these corrections: one power 
of a s comes from the Darwin interaction, Hi) ar , and the other power of a s (including the ln(Ci?a s ) 
terms) originates from the 0(a 2 ) contribution of the vacuum polarization function U Thr F " . As 
mentioned in the previous section, the latter contribution is of long-distance origin and therefore 
governed by the soft scale. In contrast to the pure Coulomb interaction, 1/p 2 , the Darwin interaction 
is a constant in momentum space and consequently sensitive to both low and high momenta. But 
based on our previous observations of the domination of long-distance effects, we can assume that 



10 It should be noted that the (1 — 4CVas/7r) correction factor of the cross section is irrelevant for shifts of the bound 
state energies because the former represents a global multiplicative short-distance factor. 

11 In pa B3] a proof for saturation is only given for the cross section above the threshold point. An analogous proof for 
highly excited stated or the cross section slightly below the threshold point does, to our knowledge, not exist in literature. 
Such a proof is, however, much more more difficult due to the breakdown of time-independent perturbation theory for 
the logarithmic kernel SV(r) ~ ln(r)/r for high radial excitations (see e.g. 123]). Nevertheless we find it plausible that 
saturation also takes place slightly below the threshold point because the cross section at the threshold point q 2 = 4Mg 
should be well denned. 
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the scale of the strong coupling in the Darwin interaction should also be the soft scale rather than 
the heavy quark mass. The size of the strong coupling governing the 0(C F a 2 ) Darwin corrections of 
eqs. fl35|) and (|37|) can therefore be estimated via the self-consistency equation 



«.s 



a s (C F a s M Q ) 



(42) 



which leads to a s = 0.13 — 0.16, 0.25 — 0.38 and 0.34 — 0.59 for the top, bottom and charm quark 
systems, respectively. The latter ranges are obtained by using the MS definition for the strong coupling, 
the one-loop QCD beta-function and a s (M z = 91.187 GeV) = 0.125 and by taking twice and half the 
argument of the strong coupling on the r.h.s. of relation (fl2|). Further, the mass values in equation (42) 
have been taken to be the pole values. For the quark (pole) masses we have chosen M% = 175 GeV, 
Mb = 5 GeV and M c = 1.7 GeV. The reader should note that the prescription given above to calculate 
the size of the strong coupling is far from being unique. Depending on the choice of the definition of 
the strong coupling, the quark mass values or the number of loops in the QCD beta-function, larger or 
smaller values for a s might result. This dependence on the prescription is particularly strong for the 
charm systemPj. As a consequence the theoretical uncertainties quoted in this work should be more 
understood as good guesses rather than strict theoretical limits. However, we think that the ranges 
of the strong coupling given above are good enough in order to illustrate the impact of the 0{C F a 2 s ) 
Darwin corrections in particular for ti production in the threshold region. We also want to emphasize 
that our conclusions for the perturbativity of the different heavy quark systems do not depend on 
different prescriptions for the strong coupling. 



n 


1 


2 


3 


4 


oo 


Ait 


-0.05/ - 


-0.04 


-0.04/ - 


-0.03 


-0.03/ - 


-0.02 


-0.03/ -0.02 


-0.02/ -0.02 


Abb 


-0.20/ - 


-0.11 


-0.09/ - 


-0.06 


-0.06/ - 


-0.05 


-0.05/ -0.04 


-0.02/ +0.01 


ACC 


-0.34/ - 


-0.17 


-0.10/ - 


-0.09 


-0.06/ - 


-0.01 


-0.05/ +0.03 


-0.01/ +0.15 



Table 2: The relative 0(C F a 2 s ) Darwin corrections to the moduli squared of the I = bound state 
wave functions A^. n are given for the tt, bb and cc system, respectively. Displayed are the smallest 
and largest values for the range of a s values given below eq. (42) for the radial quantum numbers 
n = 1, 2, 3, 4 and for n — > oo. 



In Table ^ the smallest and largest values for the relative 0(C F a 2 ) Darwin corrections to the 
moduli squared of the I = bound state wave functions A^ n = 5 | v I / S CD (0)|| )ar ./|^ CD (0)| 2 for the 
different heavy quark systems are displayed for the ground states (n = 1) and the first three radial 
excited states (n = 2,3,4), employing the ranges for the strong coupling as given below eq. ([42]). 
For illustration the corresponding value for (n — ► oo) is also presented. The absolute values of the 



1 As an example, using the two-loop QCD beta-function results in a s — 0.13 — 0.17, 0.27 — 0.44 and 0.38 — 0.76 for 
the top, bottom charm system, respectively. At this point it is clearly obvious that the situation for the charm system 
is rather hopeless as far as the question of perturbativity is concerned. 
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Figure 2: The relative 0(Cpa 2 ) Darwin corrections to the ti production cross section in the threshold 
region for the cases a s = 0.13 (solid lines) and 0.16 (dashed lines) for stable (a) and unstable (b) top 
quarks. The circles in Figure (b) indicate the location of the IS* Coulomb energy level. 



corrections to the ground states amount to 4% — 5% for the ti, 11% — 20% for bb and 17% — 34% 
for the cc system. It is an interesting fact that for the bb and cc systems the size of the corrections 
is rapidly decreasing for higher excited states. In particular, the sensitivity of the corrections to the 
different values of a s seems to be surprisingly small for the excited states in the bb and cc systems. 
We will come back to this point later. 

In Fig. ||a and ||a,b the relative 0(Cpa 2 ) Darwin corrections to the (stable) quark-antiquark 
production cross section ARqq = ^R^hr Dar/^-Thr are displayed above the threshold point for the 
three heavy quark systems in the range < (3 < Cfo s . (For ti production this corresponds to the 
energy range < E < 5 (8) GeV for a s = 0.13 (0.16)) The solid (dashed) lines correspond to the lower 
(upper) a s value given below eq. fl4"2"|). For the ti system the size of the relative corrections is quite 
stable between —1.9% and —1.0% with the tendency to decrease in magnitude for larger velocities. It 
is striking that the dependence of the corrections on the changes in the a s value is weaker for larger 
velocities (0.3% for /3 = and 0.05% for (3 = Cpa s ). For the bb system the corrections vary between 
—2% (lower value) and +5% (upper value) where the larger values occurs for larger velocities. In 
contrast to the top system the dependence of the corrections on the changes in the a s value (3% for 
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Figure 3: The relative 0(C F a 2 ) Darwin corrections to the bb (a) and cc (b) production cross section 
in the kinematical region < (3 < Cf(x s above threshold. The solid line corresponds to a s = 0.25 
(0.34) and the dashed line to a s = 0.38 (0.59) for the case of bb (cc) production. 



(5 = and 5% for (3 = CfC( s ) increases for larger velocities. This indicates that the perturbative 
approach employed in this work works better for the ti than for the bb system. For the cc system, 
on the other hand, the dependence on the changes in a s are tremendous. Depending on the size of 
the coupling the corrections vary from —1% to +15% for j3 = up to +3% to +26% for j3 = Cfo s , 
drawing a rather uncomfortable picture for the perturbativity in the charm system. For the case of ti 
production we have also plotted the corrections for a finite width Ft = 1.5 GeV (see Fig. |2]b) in the 
energy range —10 GeV < E < +10 GeV in order to demonstrate the impact of the 0(C F a 2 ) Darwin 
corrections in the presence of the large top width. This has been achieved by the naive replacement 
E — > E + iT t in eqs. (|37|) and (|38|). We want to mention that the inclusion of a finite width by this 



naive procedure does not represent a consistent treatment at the 0(oq) accuracy level. However, we 
find that this approach is justified here in order to demonstrate that the typical size of the 0(C F a 2 s ) 
Darwin corrections is not altered if the top quark width is taken into account. In this case the relative 
0(C F a 2 ) Darwin corrections amount to —6% to —2% around the IS peak and to —2% to —1% for 
higher energies. For a rigorous treatment of the corrections due to the off-shellness of the top quark 



we refer the reader to [27] and references therein. 
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Although the 0(C F a 2 ) Darwin corrections discussed above represent only a small part of the 
full 0(a 2 ) corrections, we believe that their size can be taken as an order of magnitude estimate 
for the sum of all 0(a 2 ) corrections. We therefore have to face the questions whether or how far a 
perturbative expansion in the strong coupling in the threshold regime makes sense. Because we take 
the position that one should not automatically reject the possibility of a perturbative treatment of 
long-distance effects, we think that the 0(C F a 2 ) Darwin corrections determined in this work provide 
us with important hints toward an acceptable answer to this fundamental question from the point of 
view of perturbation theory itself. There is no doubt that perturbation theory in the strong coupling 
is still viable for the ti system. It has been shown in |20| by using more general arguments that the 
large top mass and width serve as a screening device which protects the ti properties in the threshold 
region from the influence of non-perturbative effects making the ti system into the "hydrogen atom 
of the strong interaction" . Thus a perturbative treatment of the ti system should exhibit an excellent 
convergence. This is consistent with the observations from the previous discussions showing that the 
0(C F a 2 ) Darwin corrections for the top system are at the level of a few percent for the most of the 
threshold regional. This, on the other hand, allows us to conclude that the theoretical uncertainty 
of all present analyses for the total ti cross section in the threshold region is at the few percent 
level, because no full 0{a 2 ) treatment has ever been accomplished there. Further, the theoretical 
uncertainty of such a complete analysis would then be roughly one to two percent around the 15 
peak and below several per mille for higher energies. This can be estimated by taking the a s values 
presented below eq. ( f42| ) cubed (assuming that no scales lower than Cfcz s Mq are relevant for the 
corrections beyond the 0(a 2 ) accuracy level) and by observing the sensitivity of the 0(C F a 2 ) Darwin 
corrections to changes in the values of the strong coupling (see Fig. ^b). To achieve an accuracy much 
below the percent level at the IS* peak a more rigorous treatment of the scale in the strong coupling 
governing the 0(a 2 ) corrections would be needed, i.e. an 0(a^) calculation. 

As far as the bb system is concerned, the situation is worse than for the it system. It has been 



shown in a number of classical papers [30, 31, 32 1 that a proper theoretical description of the bottom 



system can only be achieved by taking into account non-perturbative corrections, which cannot be 



calculated from first principles in QCD. On the other hand, it has been demonstrated in [ 23 1 that 
a quite acceptable "parameter-free" description of the (S-wave, I = 0) bb bound states with low 
radial excitation is possible by using perturbative calculations supplemented by non-perturbative 
contributions in the form of the quark or the gluon condensates. However, the latter analyses (as far 
as corrections to the moduli squared of the wave functions at the origin and to the cross section above 
threshold are concerned) were essentially based on formulae including only the effects of the one-loop 
running of the strong coupling and the global 0(a s ) correction factor (1 — \Cfol s Jt:). The question 
whether the 0{a 2 ) perturbative corrections lead to a still converging series was not addressed explicitly. 
Equipped with the results for the 0{C F a 2 ) Darwin corrections, we are able to draw a rough picture 
concerning the latter question for the case of the moduli squared of the / = bound state wave functions 

13 A comparison of the size of the 0(C F o? s ) Darwin corrections with the 0(Cfc s ) corrections from the (1 — 4CVa s /7r) 
suppression factor is slightly misleading in this context because the latter represents a pure short-distance contribution. 
Therefore the 0(Cfoi s ) correction should not be included in a discussion on the convergence in the perturbative descrip- 
tion of long-distance corrections. However, for the convenience of the reader, the size of the large 0(Cfcis) corrections 
shall also be given. It has been shown in M, Eq, |29| in a two-loop analysis that the scale in the strong coupling of the 
0(C F a B ) suppression factor is e' 11/24 Mq in the MS scheme. This results in -4C f a s /;r = -20%, -41% and -64% 
for the top, bottom and charm systems, respectively, using the one-loop QCD beta-function, the pole mass values given 
below eq. (E|) and a s (M z = 91.187 GeV) = 0.125. 
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at the origin. For the ground state the 0(a%) corrections should be between 10% and 20% (where the 
actual sign of the corrections can only be determined by a complete 0(a 2 ) analysis) with theoretical 
uncertainties of order ±5% coming from the ignorance of the actual scale of the strong coupling and 
other corrections beyond the 0(a 2 ) level. This does not represent an overwhelming convergence, but 
it is acceptable compared to the precision of experimental measurements [ J33| ] and it indicates that an 
actual determination of all 0(a 2 ) correction would lead to a considerable improvement of the precision 
of the theoretical description. It is remarkable that the (D(Cpa 2 ) Darwin corrections seem to indicate 
that the size of the 0(a 2 ) corrections including their sensitivity to changes in the value of the strong 
coupling is much smaller for higher excited states (see Table 0). Here, however, non-perturbative 



contributions get more and more out of control f3l| , 32] and a complete 0(a 2 ) analysis is therefore 
necessary to give a trustworthy interpretation of this phenomenon. The latter remark is also true for 
the cc system. 

Finally, we also want to mention the cc system. In view of the 0(Cpa 2 ) Darwin corrections, 
we can expect the size of the complete 0{a 2 s ) corrections to the modulus squared of the ground state 
wave function at the origin to be at least at the level of 15% to 35% with theoretical errors which 
might be almost as large as the size of the 0(a 2 ) corrections themselves. (Again we can estimate 
the size of the corrections beyond the 0(a 2 ) level by taking the long-distance a s values given below 



eq. (42) cubed.) It is evident that in the case of the cc system the limits of perturbation theory are 
reached or even exceeded. Even with a complete determination of all 0(a 2 ) corrections the theoretical 
uncertainties would not decrease considerably, which is obviously a consequence of the large size of the 
strong coupling. We therefore conclude that it will be extremely difficult (if not impossible) to achieve 
a perturbation theory based theoretical description for the charm system with uncertainties lower than 
several times 10% if there is no (unforeseen) cancellation among different types of corrections. 

To conclude this section there is a remark in order: for the calculations of the 0(Cpa 2 ) Darwin 
corrections we used the renormalized Green function at zero distances, eq. (JT6J) , without any further 
explanation. This is slightly misleading because it implies that the 0{Cpa 2 ) Darwin corrections to 
wave functions and cross sections can be uniquely separated from all the other 0{Cpa 2 ) corrections. 
As far as the lna s contribution and the digamma term are concerned this is definitely true, but this 
is not the case for the constant terms. This is a consequence of the divergences which arise during 
the calculations and which have to be renormalized. The use of our renormalized zero-distance Green 
function represents one possible way to achieve this renormalization. Nevertheless we think that our 
approach is justified in order to illustrate the possible size of the complete 0(a 2 ) corrections. This 
view is also supported by the explicit results for all 0(Cpa 2 ) corrections to the / = wave functions 
at the origin and the cross section, which will be published shortly. However, we want to emphasize 
that the latter considerations do not affect the validity of the expressions for the vacuum polarization 
function presented in Sections ^ and |3[ There, all constants are correct due to proper matching to the 
well established one- and two-loop expressions n' 1 ) and n^ 2 ', eqs. (g) and (0). 
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5 Comment on Threshold Effects far from the 
Threshold Region 

In this section we want to comment on the use and the interpretation of the expression of the QED 
vacuum polarization function valid for all energies to 0(a 2 ) accuracy, eq. (|l7|). 

We have shown in Section y that the function A, which represents the resummed expression for 
diagrams with the instantaneous Coulomb exchange of two and more longitudinal polarized photons (in 



Coulomb gauge) , see eq. ( 10 ) , essentially has to be added to the one- and two- loop expressions for the 
vacuum polarization function in order to achieve 0(a 2 ) accuracy in the threshold region \/3\ < a. Far 
from the threshold regime, however, A represents contributions of order a 3 and higher and therefore 
is irrelevant. This is what we mean by using the term "valid for all energies to 0(a 2 ) accuracy". - 
But not more! 

At this point the reader might be tempted to apply formula (|17|), as it stands, for an energy 
regime far from the threshold in the believe A would represent higher-order information which should 
improve the accuracy of the one- and two-loop expressions calculated in the framework of conventional 
perturbation theory. Let us illustrate such a scenario for the energy regime where q 2 is close to zero. 
In this kinematic region, formula fll7| ) can be expanded in terms of small q 2 . Taking into account only 
the first non- vanishing contributions in q 2 jM 2 and including only contribution up to 0(a 3 ) the result 
reads 

rrOCa 2 )/ 2^ 9 2 ^0 ra\ 1 q 2 /a\2 41 q 2 3 vr 2 q 2 4 

where the numerical coefficient of the 0(a 3 ) coefficient is 7r 2 /48 = 0.21. The corresponding multi-loop 
expression including also the first non- vanishing three-loop coefficient H reads 



QED {Q ' Vj 15 M 2 + V J 162 M 2 


<)* 


8687 vr 2 /l 1 \ 22781 " 
13824 ' 3 \8 5 / ' 27648 . 



The numerical value of the constant term in the brackets is 0.32. It is evident that the 0(a 3 ) contri- 
butions which come from IIqed and therefore contain information on the formation of positronium 
bound states are much larger than the three-loop contributions. The ratio between the former 0(a 3 ) 
contributions and the three- loop result even diverges for q 2 — > 0. The overall conclusion of this sce- 
nario would be that threshold (and therefore long-distance) effects dominate not only in the threshold 
regime but also the energy region \q 2 \ <C 4M 2 . This is obviously wrong! The "threshold effects" 



in eq. ( f43|) contradict the Appelquist-Carrazone theorem [34] and even represent contributions non- 



analytic at q 2 = 0. The solution to this apparent paradox is that Hq^u only describes the vacuum 
polarization function to O(o 2 ) accuracy. All contributions of order a 3 or higher have to ignored and 
do not represent proper higher-order contributions. This means that the contributions of the function 
A are necessary to achieve 0(a 2 ) accuracy in the threshold region, but should be neglected if the 
vacuum polarization function has to evaluated far from the threshold point. 

To make the latter point more explicit, let us imagine that the analytical form of the complete 
three-loop contributions to the vacuum polarization function would be known for all energies (in the 
same sense as the they are known for the one- and two-loop contributions, IT 1 ) and IT 2 )). We then 
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could try to determine the expression for the vacuum polarization function valid to C(a 3 ) accuracy 
for all energies in the same way as we have determined IIq E d , which is valid to 0(a 2 ) accuracy for all 
energies. This would be achieved by matching the three-loop expression for the vacuum polarization 
function to the corresponding 0{o?) formula calculated in NRQED in the same way as presented in 
Section |2[ The vacuum polarization function valid to C(a 3 ) accuracy for all energies would then have 
the formP] 



+ A(a,/3) -a 



IT 2 



24/3 



+ A(a,/3). (45) 



J" - 2 



In the second line of eq. (Eq) the contribution a [ i ^-3 ] has to be subtracted in order to avoid double 
counting in the threshold regime since 

3 3 n(3)(, 2 )^< 1 a3[^] + ^). (46) 

It is therefore clear that far from threshold the second line of eq. fl45| ) only contains contributions 
of order a 4 and higher (see eq. (19)). Expanding now IIqbS for small values of q 2 would give a 
result identical to the three-loop expression, eq. (|44|). The large non- analytical (D(a ?y ) contribution 
which appeared in eq. fl43| ) would be gone. It is obvious that this large contribution originates from 
the leading non- vanishing term of IT 3 ) in an expansion for |/3| <C 1 evaluated for small q 2 . These 
contributions survive in IIqed , eq. (|i~7|), but are cancelled in IIqed , eq. (45). Using the same line of 
arguments it can easily be shown that all contributions of the function A would be cancelled if formulae 
for the vacuum polarization function with successively higher accuracy would be determined. 

The physical picture behind this cancellation can be drawn as follows: the contributions in 
function A are generated by vacuum polarization diagrams with the instantaneous Coulomb exchange 
of two and more longitudinal photons, where the latter are defined in the Coulomb gauge. In the 
threshold region the exchange of these longitudinal photons represents the dominant effect, whereas 
all the other interactions, for simplicity reasons called "transverse" in the following, can be neglected in 
a first approximation. Although this approach is obviously not gauge invariant from the point of view 
of full quantum electrodynamics, the violation of gauge invariance is vanishing in the non-relativistic 
limit. This is not true, however, far from the threshold point. There, contributions from longitudinal 
and transverse photons are equally important. Their individual sizes are extraordinarily large but 
with different signs. Therefore, adding the transverse contributions to the contributions of function 
A the greater part of the large corrections will be cancelled off, leaving results which can be obtained 
from conventional (multi-loop) perturbation theory. This remains true at any level of accuracy. From 
this picture it should be clear that neither effects from the formation of e + e~ bound states nor from 
the Coulomb rescattering, if the relative velocity of the e + e~ pair is much smaller than the speed of 
light, can ever lead to large corrections of the vacuum polarization function far from the threshold 



14 In eq. (|45|) A(a,/3) denotes the 0(o?) NRQED contributions, including the necessary subtractions in order to avoid 
double-counting. The actual form of these contributions is irrelevant here because we only want to discuss the large 
0(a 3 ) contributions in eq. (Kt3|). However, it is straightforward to see that A contains terms of order o? in the threshold 
regime, but is of order a 4 far from the threshold point. 
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region. There, the contributions of the function A represent unphysical (and gauge non-invariant) 
contributions which cannot even be used to estimate the size of the real higher-order corrections^. 

We would like to remind the reader that the previous arguments are not applicable if a high 
number of derivatives of the vacuum polarization function below the threshold region is considered, 
(d/dq 2 ) n H(q 2 ), n S> 1. In the latter case threshold effects are essential. This can be easily understood 
from the relation 

For large n and \q 2 \ <C 4M 2 the high-energy contributions in the dispersion integration are strongly 
suppressed, which leads to the domination of effects coming from the threshold region. This fact is 
the foundation of QCD sum rule calculations. At this point we would like to take the opportunity 
to comment on a recent publication where QCD sum rules have been applied to extract a s and the 
bottom quark mass from experimental data on the T resonances [||]. In this publication it is claimed 
that 0(a 2 ) corrections to the moments .M^ CD (0) have been calculated because two-loop corrections 
to the bb production cross section have been included into the analysis. It should be clear from the 
discussions of Section || that a two-loop calculation of the cross section is not sufficient to describe 
the 0(a 2 ) corrections to the cross section in the threshold region. In the analysis of || this can 
be easily seen from the fact that the removal of the two-loop contributions (after subtraction of the 
corresponding leading and next-to-leading threshold contributions) essentially has no effect on the 
results (see Table 4 in ||). The latter observation is taken as a "final test on the importance of 
higher-order corrections". However, as shown in Section |], the 0(a 2 ) corrections to the cross section 
in the threshold region are expected to be at the 10% to 20% level and will therefore have a large 
impact on QCD sum rule calculations in the large n limit. The mistake in the arguments of M is 
that it is implicitly assumed that the Sommerfeld factor, eq. (p8|), accounts for the resummation of all 
long-distance effects. Therefore all corrections to expression (|38|) should be calculable by fixed-order 
loop calculations alone. This is true for the 0(a s ) short-distance correction factor (1 — ACfol s /it), 
but this is not the case for higher-order corrections like the 0(Cpa 2 ) Darwin corrections calculated 
in Section ||. This fact will be demonstrated explicitly in a future publication, where all 0(Cpa 2 ) 
corrections to the cross section will be presented. In is it also assumed that the effects of the running 
of the strong coupling in the Sommerfeld factor can be determined by insertion of the effective running 
coupling ay, which is related to the short-distance corrections of the QCD potential p7], S3]. We would 
like to emphasize that this approach is not justified for large n QCD sum rule calculations because 
the important saturation effects are neglected in this procedure. As a consequence the calculations 
presented in j|] are not only not at the 0(a 2 ) accuracy level but also include a systematic error at 
order a s and therefore contain much larger uncertainties than presented there. The authors of || 
finally criticize an older QCD sum rule calculation by Voloshin Q on the same subject, claiming that 
in Q the magnitude of higher-order corrections was underestimated. In this point we agree with 
the authors of ||] because in Q it is assumed that 0(a 2 ) corrections have "no enhancement" in the 
large n limit and therefore should be of order 1/n. This assumption is essentially equivalent to the 
statement that all 0(a 2 ) corrections to the ground state Coulomb wave function at the origin of a 
bound bb pair should vanish. We have shown explicitly in this work that this assumption is not true 
by calculating the 0(C F a 2 ) Darwin corrections. The latter corrections amount to 10% to 20% for the 

15 If applied to QCD our conclusion is essentially equivalent to arguments employed in J33, B6| . 
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modulus squared of the ground state wave function at the origin for a bound bb pair and are far from 
being negligible. We therefore conclude that the results presented in []7| actually contain theoretical 
uncertainties at the 10% to 20% level, an order of magnitude larger than claimed there. 



6 Summary 

In this work we have used the concept of effective field theories to calculate the 0(a 2 ) corrections to the 
QED vacuum polarization function in the threshold region and to define a renormalized version of the 
zero-distance Coulomb Green function. In the framework where non-relativistic quantum mechanics 
is part of an effective low energy field theory (NRQED), long-distance effects (coming from typical 
momentum scales below the electron mass) are determined completely by classical quantum mechanics 
calculations, whereas short-distance contributions (coming from momentum scales beyond the electron 
mass) are included via the matching procedure. For the latter contributions multi-loop techniques (in 
conventional covariant perturbation theory) have to be employed. We have demonstrated that the 
effective field theory approach represents a highly efficient method to merge sophisticated multi-loop 
methods with well-known textbook quantum mechanics time-independent perturbation theory. From 
the physical point of view this is achieved because the effective field theory concept allows for a 
transparent and systematic separation of long- and short-distance physics at any level of precision. 
For our calculations we have used a "direct matching" procedure which can be applied if the multi- 
loop results to the quantity of interest are at hand. This direct matching allows for a quite sloppy 
treatment of UV divergences in the effective field theory, but is of no value if calculations of quantities 
are intended for which no multi-loop expressions are available. 

We have demonstrated the efficiency of our approach by calculating the 0(a 6 ) "vacuum po- 
larization" contributions to the positronium ground state hyperfine splitting without referring back to 
the Bethe-Salpeter equation, and by determining the 0(Cpa 2 ) (next-to-next-to- leading order) Dar- 
win corrections to heavy quark-antiquark bound state wave functions at the origin and to the heavy 
quark-antiquark production cross section in e + e~ annihilation (into a virtual photon). If the 0(C F a 2 ) 
Darwin corrections are taken as an order-of-magnitude estimate for the complete (yet unknown) 0(a 2 ) 
corrections, the typical 0(a 2 ) corrections for the tt production cross section can be expected at the 
few percent level for most of the threshold region. Around the 15 peak they might even amount to 5%. 
For the modulus squared of the ground state wave function of a bound bb pair (applicable to T(15)), 
the 0(C F a 2 ) Darwin corrections are between 10% and 20%, whereas the corresponding corrections for 
the cc system are between 15% and 35%. The uncertainties arise from the ignorance of higher-order 
corrections, in particular from the ignorance of the exact scale in the strong coupling. We conclude 
that the determination of all 0(a 2 ) corrections would represent a considerable improvement of the 
present precision of theoretical calculations to the tt and bb system in the threshold region. For the 
cc system, on the other hand, this seems to be doubtful, a consequence of the large size of the strong 
coupling. 

Finally, we have also discussed whether the formation of positronium states can lead to large 
corrections of the QED vacuum polarization function far from the threshold region and came to the 
conclusion that such corrections do not exist. 
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